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Abstract 



A general non-relativistic field theory on the plane with couplings to an arbi- 
trary number of abelian Chern-Simons gauge fields is considered. Elementary 
excitations of the system are shown to exhibit fractional and mutual statistics. 
We identify the self-dual systems for which certain classical and quantal as- 
pects of the theory can be studied in a much simplified mathematical setting. 
Then, specializing to the general self-dual system with two Chern-Simons 
gauge fields (and non- vanishing mutual statistics parameter), we present a 
systematic analysis for the static vortexlike classical solutions, with or with- 
out uniform background magnetic field. Relativistic generalizations are also 
discussed briefiy. 
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Typeset using REVT^ 



I. INTRODUCTION 



In two dimensional space, we can have particles obeying fractional statistics [|l| and in 
the field theoretic context a similar effect is generated by introducing the Chern-Simons 
(C-S) term [0 in the action. The C-S field theory was then found to be useful in describing 
the fractional quantum Hall states 0. According to some recent suggestions 0], a suitable 
generalization of it may in fact provide a unified mathematical approach to the long-distance 
physics of various quantum topological fiuids. The construction involves a set of abelian 
gauge fields (/ = 1, 2, . . . , A^) with the C-S-type term 



(here, (nij) is a real symmetric matrix) as the sole kinetic energy density for them. By 
considering matter fields with generic gauge-invariant couplings to the a^'s, we obtain a 
compound or multi-layered system which exhibits fractional statistics for the exchange of 
indistinguishable particles and mutual 'statistical' interactions between particles belonging 
to different species (or layers). Models of a similar nature have been considered also in Ref. 
^ and, as these authors emphasize, parity need not be broken in a field theory with an even 
number of C-S gauge fields. See Ref. for the application to the quantum Hall effects in 
the double-layer electron system. 

The purpose of this paper is two-fold. Firstly, we clarify the precise nature of a non- 
relativistic quantum field theory incorporating the above form of C-S-type interaction. In 
particular we find explicitly the corresponding first quantized description in the general n- 
body sector. In the latter description, fractional and mutual statistics for the particles are 
the manifestations of the Aharonov-Bohm effect involving a combination of fictitious charges 
and localized fiuxes affixed to them. This is described in Sec. II. (For a complementary 
discussion on mutual statistics from the braid group viewpoint, readers are referred to Ref. 



Secondly, in Sec. HI, we identify the corresponding self-dual system (with an arbitrary 
background magnetic field), which has a simpler mathematical structure than the generic 





i 
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case due to hidden supersymmetry [§. For instance, thanks to the supersymmetry, one 
can construct the exact many-body ground state wave function in this case. (See Ref. 
for related discussions). In this paper, however, we concentrate on the analysis of static 
soliton solutions to the classical field equations that follow, i.e., look for a generalization of 
the Jackiw-Pi solutions |TD|,|TT|. More general types of vortex solitons, which are likely to 



exhibit fractional and mutual statistics themselves, are found. The self-duality equations for 
our model share certain common elements with those for self-dual non-abelian C-S vortices 
discussed recently [l^JT^. For instance, the Toda-type equation 



VMn|0p(r)|2 = -if,,,|0,,(r)|2 (1.2) 

has a prominent role in both cases. But note that, in our model, K = {Kppi) need not be 
equal to the Cartan matrix of a certain Lie algebra and very little is known for this case. 
Concentrating on the case of two C-S gauge fields, we will present a fairly detailed study 
on the nature of possible self-dual vortex solutions, with or without uniform background 
magnetic field. For the case with non-zero magnetic field the self-duality equations were 
also touched upon in Ref. [0, but there is only a minimal overlap between the latter work 
and ours. 

Fully relativistic self-dual systems, including several abelian C-S gauge fields, are also 
possible and we briefly discuss them in Sec. IV. Nature of static soliton solutions is discussed 
in the special case of these models. Section V contains a summary and discussion of our 
work. There are two appendices. In the first we provide the derivations of certain formulas 
appearing in Sec. II. The second appendix contains the index theorem analysis for the 
self-dual system being treated in this paper. 
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II. NON-RELATIVISTIC QUANTUM FIELD THEORY 



Choosing the basis where the matrix (k/j) (see Eq. ( p. .11 )) is diagonal, let us consider 
the non-relativistic C-S gauge field theory defined by the Lagrangian density^ 
TV . Mr / a i ^ \ 



1=1 p=l I \ui n 



C/(*'.*), (2.1) 



where (p = 1, . . . , M) denote M different bosonic (fermionic) fields satisfying the equal- 
time (anti-)commutation relations 

[%ir,t),%ir',t)]^ = [vl/t(r,t),vl;t,(r',t)]^ = Q, 

[vl/p(r, t), vl/t,(r', t)]^ = V<5'(r - r). (2.2) 

[The suffix — (+) refers to the commutator (anticommutator)]. This system possesses [f/(l)]^ 
local gauge invariance in connection with independent C-S gauge fields a^, and we have 
included the external electromagnetic field A'^^ for the sake of generality. Excluding the 
pathological cases, we will below assume that all k/'s are non-zero, and M > N, i.e., the 
number of the C-S fields does not exceed that of the matter fields. Also, for definiteness, we 
shall take the potential t/(^I/^ ^I/) to have the general form 

[/ = E^p(r,t)^J(r,t)*p(r,t) 
p 

+ \Y.j d^' t)4(r', t)V,,,{v - v')m,iv\ t)m,{v, t), (2.3) 

where Vppi{v — r') = Vpip{v' — r). 

The stationary action principle for varying yields the Gauss laws (z, j = 1 or 2) 

^ e^^W.a] = -- E i'^l'^P ■ (2-4) 

Kl p 



Note that /i, z^, A = 0, 1 or 2 and our (2+l)-D metric is given as rj^u = diag.(— 1, 1, 1). 
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These, together with the Coulomb gauge conditions Vjof = 0, then determine a^(r, t) in 
terms of the matter densities Pp{r,t) = \&j,(r, t)\l/p(r, t), 



a[(r, t) = e^^V — ^ ^ I d"^' G{r - r')p,{r\ t), 



where G is the Green's function for the 2-D Laplacian, i.e., 



(2.5) 



V'G'(r-r') = (52(r-r'), 
G(r — r') = — In |r — r'|. 

271 



Now the Hamiltonian operator of the system can be identified with 



H = Jd'rlY. (r, t) ■ n,(r, t) + U{¥, ^) \ , 



where 



np(r,t) 



he 



D^p(r,t) 



(2.6) 



(2.7) 



(2.8) 



with the fields expressed in terms of and \E' through Eq. ( ^.5[ ) 

What we have in Eq.( |2.7| ) is the properly ordered Hamiltonian, and the corresponding 
operator field equations are 



ih- 



dt 



[^p{r,t),H] 



J2 g>o(r, t)^p{r, t) - ^D2^p(r, t) + Vp{r, t)^p(r, t) 
+ J d\' ^J,(r', t)Vpp.{v - r')*p'(r', t)^p(r, t) + 7^, (2.9) 



where TZ, the quantum correction from operator ordering as first discussed in Ref. |Tl[], is 
specified as 



E 



^ 2mptc'^ \j Kj J J Y47r2|r — r'p^ 
and the operators al{r,t) represent the solutions to the equations 



(2.10) 



(A^p)^(r,t)^p(r,t)-^J(r,t)(A,^p)(r,t)l =0, (2.11) 



or, more explicitly, 



^(r, t) = Y: ^e^W, / rfV G(r - r') Jp.(r', t) 



JAr,t) = ^ [(A*p)^(r,t)*p(r,t) - vl/t (r, t)(AvI/p)(r, t) 

From the (anti-)commutation relations ( p^.2| ) it follows that 

1 



(2.12) 



[vl>,(r', t), (r, t)] = e^^V ~qiG{v - r')^p(r', t) 



(2.13) 



and, in evaluating the commutator needed to derive Eq. ( |2.9|) , we have taken (following Ref. 



Ill) that 



[^p.(r',t),a[(r 



, o)\ |r'=r 







i.e., the quantity e*"'VjG'(r — r') at the coincidence limit r' 
zero. Now we denote 



(2.14) 



r has been prescribed to be 



J 1^1 



(2.15) 



and then 



D.^p{r, t) = I V,, - ^e'^Vj ^ Ppp> J dV ^(r - r')pp, (r', t) 



--epA-(r,t) vl>,(r,t). 



E^W(r,t) = E— e^^V, /dVG(r-r')J,.(r',t), 



(2.16) 
(2.17) 



showing that the parameters g^'s, k/'s enter the field equations only through the /5pp''s. From 
this we infer that theories with different g^'s and k/'s but with the same values for the Ppp'^s 
are physically equivalent. A simple physical interpretation for the jSpp'^s will be given later. 

To see clearly the physical content of the above non- relativist ic quantum field theory, 
we will now derive the equivalent first-quantized description. For the two-particle sector of 



a one-layer system (i.e., a single matter field), this has been explicitly performed in Ref. 
Tl| . The general n-hodj sector with the Hamiltonian ( |2.7| ) will be considered here^. Let 



1$) denote any Heisenberg-picture state vector with the total number of particles equal to 



n = J2p=i TT-p- Then the corresponding many-particle Schrodinger wave function is given by 



$(rS^...,rff,...,rS-\...,rK),t)= 



M 



p=l J Up] 



f2.18) 



where the vacuum |0) satisfies the conditions \l/p(r, t)|0) = for any p = 1, . . . , M. Now, 
using the field equations (|2.9|), we have 



z4$(r« ...,e...,rr,...,r£),t) = E(o 



dt 



|-l=M/i(r«,t)...M/i(rW,t) 



k=i \ n„\ 



A+B+C+D 



L=v,,,(rr,t)...v^M(ri^;,t)} 



(2.19) 



with 



1 



E^v^.(rr,t) 

k=i \Jnp\ 



j^^/^,.,VfY^VG(r?) -r')J,.(r',t)j ^p{v^\t) ■ ■ .^p{v^lt) 



(2.20) 



5 = E 



fc=l jHpl 



^After the completion of this paper we learned that C.-L. Ho and Y. Hosotani |14| previously 
considered the Schrodinger equation for n anyons (on a torus), starting from the corresponding 
C-S field theory. Ours is more explicit and also deals with more general Hamiltonian, and so we 
include this discussion for the sake of completeness. 
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2mr 



n 2 



(2.21) 



k=i \ rip] \ 



(2.22) 



k=i J rip] 



47r2 |v.(p) _ r'|2 



Pp' 



.{T',t)^p{rP,t)\...%{r^^lt) 



where vj^''^'* denotes the derivative with respect to r^^''. 



(2.23) 



The contributions designated as C and D above have the structure of the standard 



one-body and two-body interaction terms. So we may immediately rewrite them as |T5 

AM) ^) 



C 



1 



E Vpi-f, t) + ^ E E VpA-f - 4 \ m 



(1) 



' ■ ■ ■ ' "Af ' ' ' 



(2.24) 



D 



E E 

{P,fc) (p',fe')7^(P,fc) 



$(r 



(1) 



AM) 



(2.25) 



where the sum ^ Z](p,A;) Y.{p',k')^{p,k) has the effect of taking in once every different set of pairs 
(p, k) and (p', A;'), excluding the case with p' = p and k' = k. On the other hand, we show 
in Appendix A that the contributions B and A above can in fact be expressed as 



5=E 

(P,k) 



2m„ 



he 



(p' ,k')<{p,k) "-'^ 



•$(r 



(1) 



AM) 



1 5 • • • ! '■nM ' ' 



t) 



(2.26) 



^=E E 

(p.fc) (p' ,k')<(p,k) 



ih 



nipic^ 



7{p',k') 



(p",fe")<(p,fc) 
(p",fc")#(p',fc') 
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(2.27) 



where J2{p',k')<(p,k) denotes the sum over all indices (p', k') that appear on the left of r^^"* in 



the arrangement (r^^ 



(1) 



'■■■' 711 '•■■'1 1 • • • 1 '■ k '■■■' 



(P) 



r(f), . . . , \ . . . , rW). in Appendix A 



it is further shown (after a bit involved manipulations) that the contributions A, B and D 
above combine to give a surprisingly simple expression, namely, 



A + B + D= J2 



he 



~^^p^T{''^k'\t) 



E /5pp'G(r?) - rl?')) 

,{p',k')^{p,k) 



$ri^...'r£),t) 



(2.28) 



Using this result in Eq. ( |2.19| ), we then find that the appropriate many-particle Schrodinger 
equation reads 



d 



(M) 



z^-$(ri 



'^)= E 



he '•••'^nM^ he 



+ E ypi-r,t) + ^ E E ^PP'(rr - ri'? ^) ^ $(rl^ . . . , r^^), t) 

{p,k) {p,k) {pi ,k')^{p,k) 



(2.29) 



where we have defined 



A, ^ rT-(^) r-W^ 



^{p',k')^{p,k) 



(p',fc') {p",k")^{p',k') 



(2.30) 



The wave function $ should be single-valued with respect to every particle coordinate and 
satisfy the symmetry requirement appropriate to bosons or fermions. 



$(...,rg^...,4^^^...)=±$(...,r 



(p) 



fc2 ' ■ ■ ■ ' ■''fc 



!'■■■;■ 



(2.31) 



The Schrodinger equation (|2.29|) provides an equivalent first-quantized description for 
the quantum field theory defined by the ordered Hamiltonian operator ( |2.7| ). The entire 
effect of the C-S interactions now enters through the vector potentials ^(p,fc)' which can be 
related to the induced Aharonov-Bohm-type interactions between the charge-fiux composites 
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of suitable nature. In particular, the form ( |2.3CI|) implies that the Aharonov-Bohm interac- 



tion strength between the two particles carrying the labels p and p' is equal to —jSpp'. A 
straightforward interpretation of this is as follows. In view of the fact that we had abelian 
C-S gauge fields (/ = 1, . . . , A^), we associate with a type-p particle an A^-tuple of charges 

1 N 

(Qp, . . . , Qp) and also an A^-tuple of corresponding fluxes (— . . . , — ■^) in accordance with 



the Gauss laws (|2.4|) . Then, given two particles from type-p and type-p' each, the expected 



strength of the Aharonov-Bohm interaction will be Y^i qi{——) = —Ppp', in agreement with 



the observation we just made. [Actually, with the Aharonov-Casher interaction taken as 



an additional effect to the Aharonov-Bohm interaction, the flux affixed to a p-type particle 
will have to read (—2^, • • • , ~2^)- Note that, in interpreting the many-particle Schrodinger 
equation ( |2.29| ), we are in no way bound by the Gauss laws ( p.4|) .] 



We here remark that the above interpretation, based on A^ distinct U{1) charges and 
corresponding fluxes, is not the only possible. We will illustrate this phenomenon through 
a closer look at the A^ = 1 and N = 2 cases. With A^ = 1 (i.e. one C-S gauge field only) 
but arbitrary number of matter fields, we may write Pppi = and so assign flux — ^ to 
a particle of charge qp\ here, the charge-flux ratio is necessarily the same for all particle 
species. With N = 2 (i.e. two C-S gauge flelds), on the other hand, we have the formula 

= ^ + (p,p' = l,2,...,M) (2.32) 

from which we immediately derive the results 

Pppi < PppPp'p', if > 

PL' > PppPp'p', if 1^11^2 < (2.33) 



for any given p, p'. Equation ( |2.33| ) puts a restriction on the possible values of the Ppp' for 
M > 2 ( to be realizable by a U{1) x U (1) C-S fleld theory), and here the sign of matters 
also. As we described in the previous paragraph, this system may be related to that of 
composites carrying appropriate 2- vector charges and corresponding 2- vector fluxes. But, for 
the case with < 1, an alternative, in some sense simpler, interpretation is also available. 
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Specifically, we assign (scalar) charge g„ = 7= and fiux 0„ = ±( — ^= H to 

Vl'^il 2^|ki| 2^|k2| 

a type p particle, and then the Aharonov-Bohm and Aharonov-Casher interactions between 
two particles from type-p and type-p' will have the net strength 



gp0p' + = ± 



11 22 



(2.34) 



\Ki\ \K2\ 

i.e., equal to (3pp/ (see Eq. (|2.32| )) under the restriction < 0. This shows that a multi- 
component system of charge-fiux composites, with different charge-flux ratios for individual 
components, can equivalently be represented by a f/(l) x f/(l) C-S field theory. We have an 
exceptional situation if 

Pip' - PppPp'p' = (for every p, p'), (2.35) 

and for this case the charge-fiux ratios qp/(pp become independent of p. When the Pppi's 
satisfy the conditions ( |2.35| ), the equivalent U{1) x f/(l) C-S theory constructed according 
to the above correspondence is effectively reduced to a f/(l) C-S theory; i.e., from the two 
C-S gauge fields, it is only their particular linear combination that has a dynamical role. 
[Incidentally, if k,ik,2 < 0, we will always be allowed to set thanks to the 

rescaling freedom of the C-S fields. Then we may introduce new C-S fields f , f by 

«iM = ^«^+^f), «2, = ^(t;i^)-t;f), (2.36) 

and in terms of these fields the C-S Lagrangian becomes 

= «e'^^^t;(i)9,t;f . (2.37) 

The Lagrangian of this form has been considered recently by Wilczek .] It should be an 
interesting mathematical exercise to extend the above discussion to the case of general N, 
but it is not pursued in this paper. 

As is well-known, the Aharonov-Bohm interaction affects the statistical character of the 
particles involved. To see this, observe that the vector potential in Eq. ( [^.30D is locally a 
pure gauge, 

12 



AM).(r!^...,r(^;) = -VS-''=)A, 

A = ^ E E Pry'lt.n-^ ( %, %A . (2.38) 

So, by redefining the single-valued wave function $ according to 

$(rS^\...,rW t)=e--^$(r(^...,r(iJ),t), (2.39) 

we may have the gauge potentials A(p^k) disappear in the Schrodinger equation for $. But 
the gauge function A here being multi- valued in general, $ will have to be multi-valued (for 
a single- valued $) and so satisfy highly non-trivial boundary conditions that are sensitive to 
the /3py/'s. To be explicit, consider exchanging the positions of two identical particles, say 
r^^^ and r^^\ with the function $(. . . , r^J^ . . . , r^'^^ ■ ■ ■) along a certain closed path C (see Fig. 
1). Then from Eqs. ( p.31|) and ( p.39D , the resulting expression should differ from the original 



by the phase factor -^ew^f^w ^■h^'yC) ^ where 7(C) is equal to the sum of the /5pp/'s over the 
index set (p', k') associated with the r^^ ^'s in the interior of C . [Here an implicit assumption 
is that no position variable other than r^^"* and r^,^'* takes values on C\ Also note that as we 
allow a specific position variable r^^'' to be taken along a closed path C, the initial and final 
expressions of $ should differ by a phase e'^^^'-^\ with 7(C) defined as above. If all Ppp''s 
with p ^ p' vanish, the system is that of M species of fractional-statistics-obeying particles. 
In the terminology of Wilczek 0, the /?pp''s with p ^ p' are relevant for mutual statistics, 
while the /5pp's are responsible for more usual fractional statistics. We close this section 
with the remark that, because of the highly non-trivial nature of the boundary conditions 
satisfied by $, the regular gauge description based on a single-valued function $ should be 
preferred for all more explicit studies. 



III. SELF-DUAL SYSTEMS AND SOLITON SOLUTIONS 

The system described by the field theory of Sec. II is highly nontrivial, and it is extremely 
difficult to obtain any concrete information on its behavior. Naturally, one might then 
ask whether there exists a certain special choice of the potential f/(\l/^,\l/) for which the 
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mathematical treatment of the system becomes more tractable. This leads us to consider 
the so-called self-dual system, which is based on the potential 



p,p 



(3.1) 



where B^^{r) = e*-'VjA^^(r) and dp = 1 or —1 (for each p, independently). For this choice 
of the potential, the many-particle Schrodinger equation (p.29|) can be cast into the form 



(9$ 

-2 



{p,k) \ P/ 



he ' '•••-"«^ nc 



(p,k) \ '^'^pC 



since (3ppi = (3pip and we have, thanks to Eqs. (|2.30| ) and (|2l6| 



e'WS-'^U(p,,), = - E ^P'n^t^ - rg'^) • (3.3) 

{p' ,k')^(p,k) 

Notice that -ff(ist) in Eq- (|3.2| ) has the form of the non-relativistic (many-body) Pauli 
Hamiltonian on the plane, with spins CTp = +1 or —1. There is a hidden supersymmetry in 
the system ||,^ which can be exploited to find the exact many-body ground state. But, in 
this paper, we shall direct our interest to the static solutions of the corresponding classical 
field theory. That is, we consider \E'p(r, t) to be classical c-number fields and \l/j,(r,t) the 
corresponding complex conjugates \E'*(r, t). The Hamiltonian is as in Eq. ( |2.7|) (with the 
potential U given by Eq. ( p.l|) ), but the classical equations of motion do not include the 
operator ordering term in Eq. (|2.9|) , i.e., 

^^^^ = E?W(r,t)v^p(r,t) - ^D'^p(r,t) - ^^p^^'^'^ 



(3.4) 
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where the C-S gauge fields a^(r, t) are of course supposed to satisfy Eqs. ( |2.5| ) and ( p. 121) . 



As we will show, this classical field system under suitable restrictions admits a class of 
interesting, vortex-type, static solutions, carrying non-trivial characteristics endowed upon 
them by the C-S interactions. Our work generalizes Refs. |[Tl| and ||T^, where the case of 



a single matter field with = (Ref. |n|]) or 5^=^ ^ (Ref. [0) was analyzed. It is 
conceivable that the solitonlike solutions discussed here may have significant physical roles 
as regards the nature of various topological fluids within the effective field theory approach 

ii- 

To proceed, note that the choice of the potential as in Eq. allows us to write the 

static energy functional in the form 



E 



r 

jd^T.^ I^i^p + ^^P^2^p|' , (3.5) 



(A^p = (v. - ^ E q^al - ^^CpATj %, . = 1,2) 

dropping unimportant surface terms. This is an immediate consequence of the identity 

1 



za,e^^V,{%D,^p) (3.6) 



and the relation (|2.4] ). Hence any configuration satisfying the self-duality equations 

Di^p = -i(TpD2^p , (dp = +1 or -1) (3.7a) 
Kie'^V,a] = -Y,ql%^>p (3.7b) 

V 

will have the lowest possible energy, i.e., = 0. A solution of these equations should solve 
the classical field equations ( p.4|) automatically, but there is of course no guarantee that one 
can always find a non-trivial solution to Eqs. ( 3.7a|) and (p.7b| ). Suppose that a non-trivial 



solution exists. Then we may conveniently write (i.e., work in the Coulomb gauge) 

al (r) = ~^^V,U\v) , ATir) = -e^^V.V^^ir) , (3.8) 

and introduce the functions /p(r) {p = 1, . . . , M) by 
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For the functions /p(r), Eq. ( p.7a| ) now implies that 



(Vi + 2apV2)/p(r) = 0, (3.10) 

and therefore /p(r) should be restricted to an entire function of Z(o-p) = x + icTpV, viz., 
fp = /p(-2(o-p))- Clearly, the function \l/p(r) may vanish only at the zeros of fp and let these 
zeros be at (Ri^'', . . . , R-^-*)- Also, from Eq. ( p.9|) , we have 

5:g^^f/^(r) + epV^-(r) = -ap^ln[|vl/,(r)|V|/p(^(.^))r] , (3.11) 



and combining these with Eq. ( p.7b| ) then yields the equations 



In = (Tp^jT. /^PP'l^pf - ep^" [ + E 6\r - R^^^)) , (p = 1, . . . , M) (3.12) 

r=l 



P' 



he 

where we have used the definition ( |2.15D and the relation In |/pp = 47r X)r=i ^^(^ ~ ^i^'' 



Our problem is now reduced to the study of the coupled nonlinear equations ( p.l2[ ) 



satisfied by the matter densities |\I^p(r)p. Note that, in connection with the C-S interactions, 
only the parameters f3pp' enter Eq. (|3.12|) ; this is natural in view of our general observation 
made in Sec. II. While the equations of this form appear very frequently, say, in the study of 
various integrable models, there is no systematic mathematical method for constructing the 
solutions for generic {Ppp'), M > 2. Nonetheless, the solution space is expected to have a rich 
structure. With that in mind, we shall below study in some detail the nature of solutions 
allowed when there are just two independent matter fields, i.e. M = 2, and the external 
magnetic field B*^^ is zero or at most uniform. Assuming this, we now write (\I'i, \E'2) = {4>, x) 
and present Eq. ( p.l2| ) in the form 



ni 



In = ai(Ai|0r + A2|xr - eiB'^n + 47r ^ 6\r - R,) , 

r=l 

In \x\' = a^iPiM' + ^22|xl' - 62^^") + 47r £ ^^^r - R;) , (3.13) 



r=l 



where we have denoted -^Ppp' = Ppp' and = Cp. Here it is not difficult to guess that the 
system under non-zero behaves differently from that for B^^ = 0. Also, a rather different 
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behavior is expected when the condition /3ii/?22 = Pu is satisfied. [See Eq. (|2.35|) .] Hence, 
for the respective cases decided by these factors, we will separately look for the solutions to 
Eq. (|37[3| ) below. 



A. The case with 5^=^ = and A = /3ii/?22 - /3i2 / 
In this case, Eq. ( |3.13| ) may be written in the formQ 



Sn|0|O 








= -K 








[\x\') 



(3.14) 



where 



K 



( 
V 



-o-i/5ii -cri/?i2 







Kl2^ 


J 


^K2l 


K22 J 



(3.15) 



This is a Toda-type equation and we are here interested in the regular solutions with 
= J(fr\(f)\'^ < 00 and = /(i^ |xP < cxd. As we shall see below, the characters of this 

equation depend very much on the properties of the matrix K. 

The above system is integrable if K has certain specific forms. To discuss this case, we 

rescale the matter fields as |0p = |0p/a and = {a, b: positive constants) so that Eq. 

( p.l4| ) may assume the form 

\ 



Sn|0|2 ^ 






( 






, K = 








V 



-ai/?ii/a -ai(3i2/b 
-a2pi2/a -(T2I322, 



(3.16) 



This equation will be integrable |T3| if the matrix K is identified with one of the Cartan 
matrices of the classical Lie algebra. We here recall following rank-two Cartan matrices 



(3.17) 





i 2 -A 




( 2 -2] 




( 2 
























V 









^Here and henceforth, we will often omit the 5-function terms in the equation, which are significant 
only at the zeros of <^ or 
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Comparing K with these expressions, we conclude that our system is integrable if cri = (J2 
cr(= ±1) and K belongs to one of the following one-parameter families: 



K = a 


(-2 






(-2 


A 




(-2 


'1 






, a 






, or a 








. 1 


-2 J 










. 3 





(3.18) 



[Note that, in Eq. (|3.18|) , the constant h appeared above has been set to a, |a or |a. 



These special cases are known to be also relevant for non-abelian self-dual C-S systems first 
discussed in Ref. |jl2|. We here have A = /3ii/522— /5i2 > and hence, according to Eq. ( p. 331) , 



they all belong to the case with > 0. A constructive method for the corresponding 
exact solutions is described in Ref. [|13|. An interesting property of the solutions when K is 
equal to the Cartan matrix of A2 (and probably B2 and G2 as well) is that the charges Q^, 
Qy^ are quantized in such a way that the fluxes (see Eq. ( |2.4| )) 

= j dhe'^VMW, + = -{l3iiQ<t> + /5i2Qx) . 
$^ = jdh e'^Vi{qla] + qla]) = -{P2iQ^ + P22Qx) (3-19) 



may become integer multiples of 27rhc. 

For a generic matrix K, Eq. ( p.l4| ) has been studied little so far. For example, there 
is no known criterion for the existence of a regular solution to the equation. Therefore, we 
will first try to narrow down the range of the parameters (3pp' (for given (Xi, (72) in which a 
regular solution with finite charges might be available, and then give some explicit solutions 
for certain specific cases. First of all, we recall that a solution to Eqs. ( |3.7a| ) and ( |3.7b[ ), if 
exists, should have zero energy. Hence, in view of Eq. ( p^.Tj ), there will be no regular self-dual 
solution if the potential U is positive definite. The potential is, from Eq. (|3rT|) , given by 
(here, it suffices to set mi = m2 = m since the self-duality equations are independent of the 
masses) 



4m 
1 



^?7 = criAi|0|' + (fXi + a2)^: 
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+ cr2/522|xl^- 



(3.20) 



When both ki and K2 are positive, we have (3ii > 0, P22 > and /5ii/?22 — Pu > and so 
conclude that U is positive-definite if (Xi = cr2 = 1- With ki < and ^2 < 0, f/ will be 



positive-definite if ai = a2 = —1- On tlie otlier liand, witli ki > 0, and K2 < (and tlierefore 
/5ii/522 — P12 < 0), we will have positive-definite potential if all the coefficients appearing in 
the right hand side of Eq. ( |3.20| ) are positive, viz., no self-dual solution if Ku < 0, K22 < 
and (for the case of o"iO"2 = 1) -^12 < 0. The parameter range is also restricted by the fact 
that if the Laplacian of a function / is positive for large |r|, / is asymptotically increasing. 
We may apply this with / = — In |0p or — In noting that |0p and |xP should approach 
zero as |r| 00 (and hence the functions — ln|0p and — ln|xp increase indefinitely) to 
obtain a configuration with finite charges. Then, in view of Eq. (|3.14|) , this cannot be 
realized if the matrix K is strictly negative; hence, no solution if Kpp'{= —apjSpp') < for 
every p, p'. 



For further restrictions on the parameters, it is useful to note that Eq. ( |3.14| ), now 
including the (5-function terms, can be cast as 

/{ixf/n-.k-z;i};-"\ 

/MVn^^t^M!!lU,„etA-)|,r (3.21) 

Ulxf/nriik-z;!}*"/ 

If deti^ = o'iO'2{PiiP22 — (/?i2)^} is positive, these relations show that both 

{\x\VWii\^-Z'r\f'' , MyUZi\^-Zr\f'' , , , , , r 11 • 
and should be asymptotically mcreasmg; 

{\<f>\yWU\z-Zr\}'''' {\x\yuT=i\z-z:.\f'' 

but, this cannot be the case if K12 > and ^22 < or if K21 > and Ku < 0. With 
detK < 0, we encounter a similar inconsistency if K12 < and ^^22 > or if 7^21 < 
and Ku > 0. Based on this discussion, we conclude that no solution can be found in the 
following parameter ranges: 

(i) (cTi = (72 = +1, P12 < 0) or (cTi = —(72, /9i2 > 0), with ni > and K2 > 0, 

(ii) ((Ti = (72 = — 17 P12 > 0) or (cTi = — (12) P12 < 0), with Ki < and K2 < 0, 

(iii) (Ji = (72 = o"(= ±1), K12 = —<7(3i2 < 0, and at least one between Ku{= — C/Sii) and 
-^22(= —(^^22) is positive, with K1K2 < 0, 

(iv) (TiO-2 = -1, /5i2/322 > or /5i2/3ii > 0, with K1K2 < 0. 
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We have summarized our findings in Fig. 2. The shaded regions are those excluded on 
the basis of the above arguments, viz., there does not exist a non-trivial solution satisfying 
the self-duality equations if the parameters lie in these regions. Note that exactly a half of 
the parameter space is excluded. 

Let us now discuss some explicit solutions to the self-duality equations in question. A 
glance at Eq. ( p.l4|) shows that, with the ansatz = 7|0(i")P (7- a positive constant), 



it can be reduced to a single equation for |0p as long as 7 is chosen suitably. Indeed, as we 
make the choice 

7 = ' ' , (3.22) 

J\22 — J^12 



what follows from Eq. ( p.l4| ) is just the Liouville equation for 



|2. 



VMn|0p = -^^|0r. (3.23) 



To have non-trivial solutions, the coefficient of |0p on the right hand side of Eq. ( p.23| ) 



should be negative as well as 7 > 0. For these Liouville- type solutions, the fluxes are 
quantized, i.e., = = (integer multiples of 2t[Tic). We here note that, due to the 
restrictions mentioned above, the parameter range allowing these Liouville-type solutions 
does not cover the entire unshaded region in Fig. 2. One may then ask 

(i) When the Liouville-type solutions exist, are there additional solutions distinct from 
the Liouville-type? If so, should they have quantized flux values always? 

(ii) When the parameters are such that no Liouville-type solution exists, will there be 
some solutions to Eq. ( p.l4| ) after all? 



These issues are discussed below. 

On the question (i), the existence of more general type solutions is confirmed by the 
index theorem, and our numerical study further shows that there are also solutions with 
nonquantized fiux values. The result of the index theorem analysis (see Appendix B for 
details) is as follows: the number of free parameters in the general solution, with the field 
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(x) having vorticity rii (722) and the asymptotic behavior ~ (~ is equal to 

2(ni + di) + 2(?22 + tt2), where a denotes the largest integer less than a. Here the 'vorticity' 
is equal to the total number of zeros in a specific matter field. By applying this to the case 
when the Liouville-type solutions appropriate to the values rii = n2 = n and ai = 02 = n + 2 
are allowed, we immediately see that the corresponding general solution should contain 
8(?2 + 1) free parameters. But this is precisely twice the number of free parameters entering 
the Liouville-type solutions A plausible conjecture is that, in general solution, we no 
longer have the restriction (satisfied by all Liouville-type solutions) that the zeros of are 
also the zeros of x- Furthermore, we found numerically that there exist also solutions with 
the asymptotic behaviors not given by integral power falloff, i.e., ai and 02 above need 
not be integers but vary continuously. For an example, see the next-to-next paragraph. In 
view of the relations = ±2TThc{ni + ai) and = ±2TThc{n2 + 02), these solutions will 
then have nonquantized fiuxe^. [In this regard, see also the comment immediately after Eq. 

(H)]- 

In the parameter range where no Liouville-type solution exists, analytical means are 
not available at present and we resorted to numerical analysis (assuming the rotationally 
symmetric form). For some values of parameters we succeeded in finding solutions while, 
for other values, no acceptable solution could be found. It seems that solutions exist in a 
large portion of this parameter range also, but the precise criterion on the parameters to 
have solutions is not clear yet. We also note that if a particular non-Liouville-type solution, 

^Recently some authors argued that the flux (or charge) for the solutions of the Liouville 
equation is quantized because of the inversion symmetry. This is misleading, however. Rather, 
we might as well say that the inversion transformation generates non-singular solutions because 
the charges happens to be quantized (and, correspondingly, the field has integral power falloff 
asymptotically). No useful information is gained by considering the inversion symmetry in our 
case. 
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conformal transformation [121, i.e. 



say I , is known, a family of new solutions may be obtained by considering its 

for a polynomial function 

w{z)m^{z),^*{z*))\) 

i^{z). 

Now, as a specific example, we will discuss solutions in the self-dual system with (3ii = 
(322 = but nonzero (3i2 (i.e. keep only mutual statistical interaction). This can be realized 
by the choice 



-K2 = K, 



1 2_ 1 



and then j3i2 = The self-duality equations now read 



he K 



(3.24) 



(3.25) 



choosing o"i = = — 1 so that solutions may exist. Note that if the theory is rewritten 
in terms of the fields f^,^-* and f^^^^ defined by Eq. (|2.36| ), this system is described by the 
Lagrangian density 



2mi he'^ 



d___ 
di h 
2 h^ 



I (1) 



d 



(2) 



X 



2m9 



V - — g V 

he 



X 



+ 



h h 
+ 



qq 



?\x? 



(3.26) 



2miC 2m2e) k 

The equations (|3.25| ) clearly admit the Liouville-type solutions, which are based on the ansatz 
\4>{r)\ = |x(r)|. In addition, we have found some non-Liouville-type solutions numerically, 
assuming that the fields (f){r) and x(r) are functions of r only. These are shown in Fig. 
3. In Fig. 3(a), a plot is given for a solution corresponding to rii = n2 = (i.e. zero 



vorticity for both (p and x)? the asymptotic behavior of which is determined as 



-1.34 



and Ixl 



,-3.95 



. Evidently, the fluxes are not quantized for this solution. Another plot now 



for a solution corresponding to ni = and ^2 = 1 is given in Fig. 3(b). 
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Finally note that the charges {Q^, Qx) of a vortex soliton with given fluxes 
determined by Eq. ( p.l9| ), and they will be in general non-zero. This prompts one to make a 
conjecture that the vortices described above experience fractional and also mutual statistical 
interactions analogous to those experienced by the elementary quanta in the theory. The 
most direct approach to settle this question is to derive the effective Lagrangian which 



is relevant to the slow dynamics of these vortices p^. But, in our case, this program 
is nontrivial especially because a multi-vortex type solution cannot be interpreted in an 
unambiguous way as representing a collection of some elementary vortices here. So it remains 
as an open problem. 

B. The case with B""" = and A = ^ii/322 - /Sfa = 

If /5ii/322 is equal to /3f2, two C-S gauge flelds are no longer independent and hence we are 
led to a self-dual system with two matter flelds coupled to a single C-S gauge fleld. In fact, 
we may now write A,0 = (Vi - Ax = (V» - ^g2ai)x, and pij = = 1, 2). 

Here we flnd it convenient to write Eq. ( |3.13|) using the couplings qi = ^^^g^ rather than 
Pij, viz., 

v'in\^\' = -aMqM' + q2\x\'), 

In IxP = -^2g2(gi|0|' + q2\x\') , (3.27) 

where we have chosen k < (with no loss of generality). Evidently, if qi and q2 are of 
the same sign (i.e. gig2 > 0), these equations will admit a bounded solution only for 
o"! = (T2 = +1. For qiq2 < 0, on the other hand, one can have a solution only with 
cri(T2 = — 1. To show this, suppose that ai = cr2 = 1. Then, from Eq. ( 3.27] ) (with the 
suppressed (5-function terms put in), we have 

V2 In - In = 47r g 52(r - R,) - 47r pi g 52^r - R^) , (3.28) 



and hence 
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91 792 



C- 



X 



n: 



"2 / 
r=ll 



{C: a complex number). 



(3.29) 



But, with qxqi < 0, this is impossible: the left hand side of Eq. ( p. 291) vanishes asymptoti- 
cally, while its left hand side clearly does not. The situation is analogous for o"i = o"2 = — 1. 
Hence, the choice (Jicr2 = —1 is appropriate in the case of qiq2 < 0. 

Incidentally, the above consideration is sufficient to show that there is no solution with 



finite charge if the matrix K, defined in Eq. ( p.l6| ), is equal to the 2x2 affine Cartan 





( 2 


-2^ 




( 2 


"'1 


matrix 






or 



















. [This corresponds to the case for which Eq. (|3.16| ) can 



be reduced to the sinh-Gordon or Bullough-Dodd equation, as noted in Ref. |T2[.] Indeed, 



for these particular forms for K, we have 0"icr2 > and then, by the above consideration, a 
bounded solution may be possible only with cri = (72 = +1. So choose ai = a2 = +1 and 
then we find the matrix (/3pp') described by the one-parameter family of singular matrices, 

. When translated into the above notation, this matrix form 

leads to qiq2 < (i.e. the wrong sign); hence, there should not be any bounded solution. 

Let us now study soliton solutions for some representative cases. First we consider the 
system with qi = q2 = q (and (Xi = cr2 = +1, of course), so that Eq. ( p.27|) may assume the 
form 





(-2 


2] 




(-2 


A 


a 






or a 








y 2 








1 

"2/ 



VMn|0|2 = -g2(|0r+|xr) 



V'\n\x\' = ~q\\^\^ + \X\ 



(3.30) 



The corresponding system has additional global SU{2) symmetry and may be viewed as a 
nonrelativistic version of a special self-dual model considered by Kim |^ . By exploiting this 
global SU{2) symmetry, a series of exact solutions, which are different from the Liouville- 



type solutions (obtained under the ansatz |x| 



71 



can be obtained. Specifically, we 



found that the coupled equations in Eq. (|3.30| ) have also the solutions of the following type 



(A 



12iPiz)Q'{z) - Q{z)P'{z)) 
|g|(|P(^)|2 + |Q(^)|2)3/2 



P{z) 

Q{z) 



(3.31) 
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where P{z) and Q{z) are arbitrary polynomials of z under the restriction that these two 
functions share no common zero. Note that this solution in Eq. ( p.31| ) do not satisfy the 
Liouville equation and include no rotationally symmetric configuration since P{z) and Q{z) 
have no common zero. The fluxes and for these solutions are quantized as in the 
Liouville-type solutions; in detail, for the solution ( |3.31| ), we have $0 = $^ = 27ihc{3np + 
Suq) if P{z) {Q{z)) is an np (nQ)-th order polynomial. [Since the vorticities of and x 
are equal to ni = 2np + uq — 1 and n2 = np + 2nQ — 1, this may also be written as 
$0 = <l>;^ = 2Tihc{ni + ^2 + 2)]. However, using the index theorem argument (see Appendix 
B), we know that there must be solutions other than these two types. We are also not sure 
whether or not the fluxes are necessarily quantized for all bounded solutions to Eq. ( p. 301) . 



As another case, we choose q2/qi = —2 and ai = —02 = +1. Then, according to the 
same procedure which led to Eq. (|3.29| ), we have 

^ = C*^^ , (C: a complex number) (3.32) 

and inserting this into Eq. ( ^.271 ) yields a single equation for |0|: 

In |0p = -q! I |0p - 2|Cp g7^ \l - 5[ • (3.33) 



If we here restrict ourselves to the special case 

^ = , (3.34) 

Eq. (|3.33|) is simplified as 

VHn\<P\' = -ql\^ni-2\Cm'). (3.35) 
This is identical to the equation encountered in the relativistic self-dual C-S Higgs system 



of Refs. For the latter system, there are now rigorous existence proofs for both 

topological and nontopological soliton solutions. Only the nontopological ones are relevant 
in our case, for the other class leads to infinite charge. In particular, a rotationally symmetric 
nontopological soliton solution has the behaviors 
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(for some nonnesrative integer n), as r — > 
^ (with a > n + 2), as ^ oo 

and for this solutions we find the (nonquantized) fluxes $0 = 27ihc{n + a) and = 
—2'Khc{2n + 2a). This is another evidence for our assertion that quantized flux values are 
not to be expected generally. 



C. The case with ^^'^ / 



As a uniform external magnetic field is turned on, a spontaneously broken vacuum be- 
comes possible and correspondingly we might then have nontrivial solutions to Eq. ( p.l3|) in 
the form of topological solitons. With a single matter field, an analogous phenomenon was 
noticed in Ref. [0. Assuming A = /?ii/322 — P12 (^he A = case is considered later), 
the asymptotic values of |0(r)| and |x(r)| for a topological soliton solution should be equal 
to 



= V 5^(^2261 - ^262)/ A 



= V5'"(-A2ei + Aie2)/A 



(3.37) 



Using these vacuum values, we may now rewrite the self-duality equations as (cf. Eq. (|3.14| )) 



Sn|0|2 ^ 




/ 


\<l>\'-vl^ 




= -K 










\ 





(3.38) 



with the same matrix K as in Eq. ( |3.15| ). In case a topological soliton solution is allowed, 
it will be subject to the topological quantization conditions of the form 



e''Vi{qla] + qla]) + eiB' 



e''Vi{qla] + qX) + e2B' 



2'Khcni , 
27ihcn2 , 



(3.39) 



where rii and 77-2 are integers. 
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In view of Eq. ( p.37|) , one may hope to find a topological soliton solution only when 
the parameters of the theory satisfy certain restrictions; namely, for A > (and hence 
K,iK2 > 0), one must have 



5^^(^2261 - ^1262) > , fi'=^(-A2ei + Ale2) > . 



(3.40) 



while, for A < (and hence k,ik,2 < 0), the inequality signs in Eq. ( p.40|) should be reversed. 
Aside from these, there must be some conditions involving the elments of the matrix K 
mainly. We can obtain such conditions by studying Eq. (|3.38|) in the asymptotic region. 
For this asymptotic analysis, we set 



/*(r) = \{vl - |0(r)p) , /,(r) = " I^WH 



(3.41) 



and may instead study the hnearized form of Eq. ( p.38|) , i.e.. 





U 




f<f> 


( 






= -L 








[fx) 




[fx) 


V 



VIK21 VIK22 j 



(3.42) 



Now suppose that the 2x2 matrix L can be diagonalized, i.e., SLS ^ 
some non-singular matrix S. Then Eq. ( |3.42| ) can be cast as 



Ai 
A2 



for 



\fxj 



A2 



/ ? \ / . \ 

= s 



\fxj 



\fxj 



\fxj 



(3.43) 



Here the eigenvalues Ai and A2, which have crucial importance in determining the asymptotic 
behaviors, are the roots of the secular equation 



= det(L - XI) 



A^ - {vlKn + vlK22)X + vlv^det K) . 



(3.44) 



Now take the case of K1K2 > (and so A > 0). In this case, it is easy to show that Eq. 
( p.44| ) has necessarily two real roots. Then note that, for an acceptable soliton solution, 
and above should approach zero asymptotically in such a way that the resulting soliton 
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may have finite energy. Tliis requires both roots to be negative (i.e. Ai < and A2 < 0), 
and, in view of Eq. ( |3.44| ), this translates into the following conditions: 

vlKii + vlK22 = vlaipii + vla2p22 < , 

det K = (71(72 A > . (3.45) 
Based on these, we find that the appropriate choice for a nontrivial soliton solution is 

(Ti = (72 = —1, if Ki > and ^2 > 

ai = (72 = +1, if fti < and ^2 < . (3.46) 

A similar analysis may be repeated for the case of kik,2 < 0. For the latter case, however, 
the roots of Eq. (|3.44| ) are not always real and this introduces a certain uncertain feature 



in the analysis. Nevertheless, for K1K2 < 0, we can make the following definite statement:no 
solution exists with (7i(T2 = +1- 

In addition to the above, certain (plausible) conditions can also be derived on the basis 
of the conjecture that a regular soliton solution, if exists, is likely to fulfill the inequalities 

\<P{r)\-v^<0, \x{r)\-v^<0 (3.47) 

at least for sufficiently large r. Then, accepting this behavior, it is possible to apply the 
same reasoning as in the case of B'^^ = (see the subsection A). For instance, Eq. (p. 38 



will be inconsistent with the assumed asymptotic behavior of \4>{r) \ and |x(r)| if the matrix 



K is strictly positive; hence, no solution to Eq. (|3.38|) if Kppt = CTp/Spp' > for every p, p'. 
Also, by proceeding as in Eq. ( p.21| ), we expect that a nontrivial solution may exist only 
under the conditions 

{det K){niK22 - n2Ki2} < 0, 
(det K){-niK2i + ^2/^11} < , (3.48) 

where ni (^2) denotes the vorticity of the field (f) (x)- 
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We do not know of any analytic method developed to study the system in Eq. ( p.38|) , even 
for some special K. We thus looked for numerical solutions to Eq. ( |3.38| ), while assuming 
the rotationally symmetric field configurations. Here it should suffices to say that, at least 
for certain choices of parameters (and vorticities) which satisfy the conditions given above, 
we did confirm the existence of regular topological soliton solutions. Note that, in view of 
the index theorem (see Appendix B), the existence of a particular solution with vorticities 
rii and n2 actually imply the existence of a 2(ni + 712) parameter family of soliton solutions. 

If A happens to vanish, we again have a self-dual system with two matter fields coupled 
to a single C-S gauge field. Using the same notation as in the subsection B, it will be possible 
to rewrite Eq. ( |3.38|) (with the choice k < 0) as 



V' In 101' = -(Ti 

72 i„ |,,|2 



V'lnlxl 



-(72 



qi{qM? + q2\x?)-eiB' 



(3.49) 



A precondition to have a topological soliton solution is the existence of a non-trivial (uniform) 
vacuum solution, and in the present case this will be true only when the coupling parameters 
satisfy the relation 



ei 

qi 



62 
?2 



(3.50) 



Here a particularly simple case is obtained for ei = 62 = e and then, thanks to Eq. ( p. 501) , 
qi = q2 = q- For this special case, the system has in fact a global SU (2) symmetry and this 
is also manifest in the self-duality equations 



V' In 101' = -ai 
72 i„ I, ,12 



V'lnlxl 



-0-2 



+ \x\')-eB^ 



(3.51) 



This system admits a topological soliton solution if we choose ai = (T2 = —1 and eB^^ > 0. 
Then, Eq. ( |3.51D becomes identical to the self-duality equations found in the relativistic self- 
dual Ginzburg-Landau model with the so-called semilocal symmetry . General solutions 
for this model are discussed in those papers to which readers are referred. 
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IV. RELATIVISTIC GENERALIZATION 



Here we will introduce the relativistic self-dual U{1) x f/(l) C-S system (with no external 
magnetic field, for simplicity) and then study the static soliton solutions in the model. If one 
wishes, one may view this investigation as a direct generalization of the model considered 
in Ref. ||2^. Our model contains two complex scalar fields and x? and is described by the 
Lagrangian density 



where 



(4.1) 



D^X = [^M - + qlal)]x 



(4.2) 



and f/(</), x) is to be chosen shortly. [In this section, we set c = ^ = 1.] The Gauss laws read 

1 



1 

/t2 



(4.3) 



with the currents = -i [0*^)^0 - (L>^0)*0] and = -i\x*D^x - (^^x)*x]- 
The static energy functional is given by 



4|xl 



1 2 1 2 



1 2 1 2 



(4.4) 



where, on the second line, we have used the following relations (derived from the Gauss laws 
([4.3]), assuming time- independent fields): 

iqim)b^ - {qln^W 



qlal + qlal 



2{qW2-qkl)\<P? 
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(4.5) 



Now suppose that the potential U has the form 



U 



{a,f3 



(4.6) 



(o r (a ) 

where ap = +1 or —1 {p = 1, 2), the Ppp'^s are defined as in Eq. ( p. 151) (i.e., Pu = -^7- + ^^, 

11 2 2 

P12 = + etc.), and c, c' are arbitrary constants. Then, with the help of the identities 
analogous to Eq. ( p.6|) , it is straightforward to show that the energy functional in Eq. ( f4.4|) 
can be rewritten as 



E 



\ qm - Qki J 



12 12 ~'~ 



+ \Di(j) + iaiD2(p\^ + \Dix + «c^2^2Xp 



{^i/92i(|0r-c2)+a2/322(|xr-c'')} 



(4.7) 



Hence, for the theory defined by the Lagrangian density ( [4.1|) with the potential ([4.6|) , there 
exists a Bogomol'nyi-type bound |28| for the static energy (which is nonnegative) , 



(4.8) 



where the fiuxes and $^ are defined as in Eq. ( p.l9|) . Since E must be nonnegative, this 
bound is meaningful only for a positive value of cric^$0 + a2c''^^x. 

Of particular interest here are the solutions saturating the above Bogomol'nyi bound, 
which are realized if and only if all squared expressions appearing in the integrand of Eq. 
( [4.7|) vanish identically. This gives rise to the following self-duality equations: 



Di(f) + i(TiD20 = , Dix + icr2D2X = 



(4.9) 



and 



^ vlU,|2r„ o (\A\'2 : „ o /u i2 



(12\X\ [f^l/?21 

K2 



■ c') + a2f322{\x\' - n] 
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--ql\x\'WiP2im' - c') + a^Mh? - c")] . (4.10) 
We expect that, for the parameters in some range at least, this system of equations admit 



both topological and non-topological soliton solutions just as in the model of Refs. |l24| , |25 
To analyze these equations, we may again write (r) as in Eq. ( p.8|) and then, for the 
functions /i(r) and /2(r) defined by 

0(r) = e— iW^'+9fc^n/^(r) , ;^(r) = e-'^2^'?2^'+52^'>/i(r) , (4.11) 

Eq. ( |4.9|) reduces to the statement of complex analyticity, i.e., /i = fi{z(„^)) and /2 = 
/2(-2((T2)) "with Z(o-p) = X + icTpy. At the same time, we use Eq. ( |4.11| ) to express 6^(r) and 
6^(r) in terms of |0| and and then combine them with Eq. ( |4.1CI| ). The results are the 



following equations^: 



|2 „/2xj 

ni 



V^n = 4ai {/3n|0PMii(|0r - c') + ^2^2(1x1' - c' 

+A2|xPM2i(|0P - c') + cy^Mh? - c'')]} + 47r ^ ^^(r - 

r=l 

In = 4^2 {/52i|0PMii(|0p - c^) + o^2/3i2(|xP - c'')] 

+/322|xn^i/52i(|0P - c^) + 0-2/322(1x1' ' c'')]} + 47r ^ ^^(r - R;) . 

r=l 

(4.12) 

We have here assumed that the fields (j) and x have zeros at (Ri, . . . , RnJ and (R'^, . . . , R^^), 
respectively. Notice that there is a certain similarity between Eq. ( [4 .121 ) and Eq. ( p.l3| ). This 



is not surprising since one can recover the model discussed in Sec. Ill as the non-relativistic 
limit ill]] of this relativistic theory (restricted to the nontopological soliton sector). 



^Notice that the q's, ni and H2 enter Eq. ( [4.12| ) only through the quantities /3pp'. This is an 
expected result even in the present relativistic case, for the given equations should remain the 
same under suitable hnear transformations on the C-S fields. 
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A general investigation on possible solutions to Eq. ( fl.l2| ) is beyond the scope of this 
paper. We will below concentrate on a particularly interesting special case, the self-dual 
system with Pn = (322 = and ^ 0. Note that we studied the non-relativistic model 
under the same condition in Eqs. ( |3.24| )-( |3.26| ). Choosing the parameters as in Eq. ( p.24[ ), 
Eq. ( f4.12| ) assumes a much simpler form, viz.. 



VHn\<f>\' = ^\xm{'-c') 



VHn\x\' = ^m\x\' - c") , (4.13) 

with the 5-function terms not written out explicitly. The Lagrangian density for this system 
reads 

'''"l0r(lxP-c'^)^-^lxP(l0p-cY (4.14) 



with the Gauss laws (for time-independent fields) given by 



,/2 



e^^V.^;f = -^J° = ^t;«°|0P- (4-15) 



There are two distinct classes of solutions to the self-duality equations. The first is a topo- 
logical soliton with the asymptotic behavior 

^oo: MM^i, m^i, (4.16) 
c c' 

and the fluxes are quantized for this solution, i.e., $0 = g J(Pr e^^'ViVj^^ = 27ml and $^ = 
q' Jd"^ e'^^ViVj^^ = 27^2 (^i, integers). The second class is a nontopological soliton with 
the asymptotic behavior 

, , , ,, const , , const /. 

I'^WI-T^, IxWI-^, (4.17) 

(«! and a2 are real numbers larger than 1). For this nontopological soliton the fluxes are 
not quantized: we here have the formulas $0 = 27r(rai + ai) and $^ = 27r(?T,2 + 0^2), when 
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the field has vorticity ^1(722). A topological soliton with = 27rni and = 27rn2 
may conveniently be visualized as an assembly of |r;,i| '0-vortices' with respective centers at 
the zeros of (p and 1 vortices' with respective centers at the zeros of x- 

Notice that, for the above soliton configurations, the charges = J(frJ^ and = 
J (Pr J° are simply related to the fluxes as 

$^ = -— Qx, $x = -— Q^, (4.18) 

due to the Gauss laws (|4.15|) . This relationship suggests the existence of mutual statistical 
interaction between 0- vortices and vortices; but, an assembly of 0- vortices (or vortices) 
only will show no peculiar statistical effect. This conclusion is further supported by calcu- 
lating the angular momentum J = Jd"^ e'^^XiT^\ where T^^ denotes the momentum density 
in the theory. In fact, at least for a spherically symmetric solution based on the form 

0(r) = /(r)e^"^^ ^(r) = ^7(r)e^'^^^ 

) (r) = e'^ 4 [^^'^ (r) - ^1] , q'^?^ (r) = e'^ ^ [h^^^ (r) - ^2] (4. 19) 



with 



h^^\oo) = h^'^\oo) = , (topological soliton) 

(nontopological soliton) (4.20) 

a simple calculation gives the result 

TTK / 1 • IN 

nin2 , (topological) 

J=\ (4.21) 

{nin2 — 0102) , (nontopological) . 

qq' 

Thus, individual 0- or x-vortices do not carry angular momentum. On the other hand, a 
composite of a 0- vortex and a x-vortex each has a nonvanishing J, and this is an anticipated 
result in a system with mutual statistical interaction. 

The general solution to the given self-duality equations is difficult to obtain, although 
certain subclasses of solutions can be readily identified in terms of those of the previously 
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known self-dual system. By evaluating the index of the differential operator associated with 
the appropriate fluctuation equation, the number of free parameters entering the general 
solution with given values of $0 and is determined as 

{2ni + 2^2 , (topological) 
(4.22) 
2ni + 2ai + 2^2 + 2^2 , (nontopological) . 

Here the general topological soliton solution with ^2 = but rii ^ (or, if one wishes, 
with 77-1 = but n2 7^ 0) is easy to describe — one may set x(r) = c' and, in view of Eq. 
( [4.13| ), just choose 0(r) to be a solution to the familiar equation from the study of the 
Ginzburg-Landau-type model E^, 



In 101^ = 1^^(101^ -c^). (4.23) 



In this case, it follows from Eq. (|4.15| ) that v^^^^{r) = Vj^\r) = 0, while 



^^%^, e'^^ViVj (r) 7^ 0. Another subclass of topological or nontopological soliton solutions 



are obtained by setting |x(i")P = (7) |0(r)P, and this of course corresponds to the case 
with rii = 77-2 (and ai = 02). For the latter, the two equations in Eq. (|4.13|) collapse to one, 
namely. 



v^ini0r = ^^^i0r(i0r-c^), (4.24) 



the form of which matches precisely the corresponding equation encountered in the study 
of the 'minimal' self-dual C-S Higgs model |2^. But this does not comprise the full general 
solution in the sector specified by rii = n2 (and ai = 02). The number of free parameters 
which enter the solution based on Eq. ( [4.24]) (as calculated in Ref. |^) is just a half of 



the value given in Eq. (|4.22|) . This may be understood by observing that the basic unit in 



a solution satisfying the condition |xp = (^^^ |0p is assumed by "a 0-vortex on top of a 
X- vortex" , while the index theorem suggests the existence of more general solution in which 
0-vortices and x-vortices serve as separate units. 
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V. SUMMARY AND DISCUSSION 



The precise nature of the Schrodinger quantum field theory with general [[/(l)]^ C-S 
interations has been clarified, a novel feature being the existence of mutual statistical in- 
teractions between distinguishable particles. Then, for the corresponding self-dual models 
with two matter fields, we investigated the structure of classical soliton-type solutions to the 
static equations of motion, with or without uniform external magnetic field. In particular, 
to obtain a system which admits non-trivial soliton solutions satisfying the self-duality equa- 
tions, we derived a set of necessary conditions for the parameters of the theory. While our 
self-duality equations reduce to the Toda-type equations or their generalizations, the matrix 
K in the equations is not necessarily equal to the Cartan matrix of a certain Lie algebra. For 
some special cases we exhibited soliton solutions in a more explicit way. Soliton solutions 
in a relativistic self-dual system with two C-S gauge fields were also discussed briefiy. We 
conjectured that these solitons exhibit mutual (as well as fractional) statistics. 

Some comments are in order. First of all, it is intriguing that the Toda equation retains 
some of its interesting mathematical properties (e.g. the existence of multi-soliton solutions) 
even if its structure is suitably altered. Aside from that the matrix K in Eq. ( p..2|) need 
not have a group theoretical origin, we saw this phenomenon realized when we add constant 
terms on the right hand side (as in Eq. ( |3.38| )) and also quadratic terms in the densities 
(as in Eq. (|4.10| )). Quite possibly, certain universal mathematical structures might exist 
behind all these models. Also desirable will be to clarify further various physical properties 
(e.g. statistics) of the vortex solitons discussed in this paper and to study their possible 
roles in the real physics of multi-layered Hall media. Another fruitful line of research is 
the quantum mechanical investigation of our model Hamiltonian in Eq. (|3.2|) . We noted 
already that, by exploiting the supersymmetry in this system, it should be possible to find 
the corresponding exact ground state and also their degeneracy. Just as in the case of a one 
layer system 0, this investigation might yield some valuable insights in understanding the 
multi-layered fractional quantum Hall effects. 
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APPENDIX A: 



Here we will first explain how Eqs. (|2.26|) and ( p. 27 ) are derived, and then go on to 
establish the expression ( 2.28|) . For the contribution B defined by Eq. ( |2.21 ), one may 
repeatedly use the relations hke 



2m. 



2mr. 



he 

* he 



p' 



he 



(p) 



(Al) 



which follows from the noncommutativity of pp'{r',t) and \E'p(r^^]^, t). Once all the field 
operators on the left of squared differential operator are relocated to its right by this proce- 
dure, one readily recognizes the expression in Eq. ( |2.26|) as a consequence of the definition 
( p.l8| ) and the fact that (0|pp'(r', t) = 0. For the contribution A, more steps are necessary to 
derive the result ( |2.27| ) from the form in Eq. ( |2.2CI| ). Here, using the commutation relations 
]2D and ( |2.13|) , we first observe that 
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ih 



rripC 



T)' y 



(A2) 



where D^^'''^ is defined as in the case of V-'^'^^ If we further let the operator \I/p(r^^2i ^ct 
on the expression (|A2|) from the left, the result can then be written as 



ih 



rUpC 
ih 



TTlpC 



TUpC 



J 



ih 



nipC 



(A3) 



-Y.^2p>''^f'^ I d'r'Gir^^ - r')6^'Vr^[G(r?), - r') + G{r^l, - r')] 

p' p 

^l,{r',t)^p,iv',t) 

+ ^^e^^Vf^) [d\'G{v^^-v')JpH{r',t)] ^p{r,.,,t)%{/^l„t), 
p' ^ / 

where we have again used Eq. (|2.13| ). Now it is not difficult to infer that as analogous steps 
are repeatedly all the way, the final result should be the expression ( |2.27|) . 

To show that the sum of A, B and D can be expressed as in Eq. ( |2.28|) , we proceed as 
follows. We begin with the trivial observation, ^ = E "I" E ' '^^^ 

{p' ,k')yt{p,k) {p' ,k')<{p,k) {p',k')>{p,k) 

the right hand side of Eq. ( p. 28 ) into the form 



ih 



{p,k) ^P^ 



E Ppp'G(jc^k^ -^t'^^ 



\{p' ,k')>{p,k) 

^'f^'^ - y/^^^"^ ( E ^P"G{rt' - ^'') 1 - i-^P^n-T^t) 



Sp) Jp"h 



y{p",k")<{p,k) 



I 

he 



Sp) 



■*(rS^...,C\t) 



(M) 



+ E E E ^.PPP'PPP" - r^))) 

(p,fc) ip',k')>(p,k) ip" ,k")>{p,k) ■^'"'P'- 
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(A4) 



The second term in this expression can then be rewritten as 



{p,k) {p',k')<{p,k) ^P"^ 



(1) 



he 

AM) 



y(p" ,k")<{p' ,k') 



(A5) 



while the last term is equal to 

y y , 

^ ^ 2m ic^ 

{p,k) {p' ,k')<{p,k) ^" V- 



,(1) 



EE E (e-vf''G{vt' - r^P 

(p,k) {p\k')<{p,k) (p",k")>(p',k') "V- 
(p",fe")>(p,fc) 



(A6) 



The first term in Eq. ([A6|) evidently coincides with the contribution D shown in Eq. ( p. 251) . 
On the other hand, the second term in Eq. (|A6|) and the expression (|A5|) combine to yield 
the contribution A in Eq. ( |2.27|) . The equation ( p.28| ) has been established now. 



APPENDIX B: 



As regards the solutions to the self-duality equations ( |3.7a| ) and ( p.7b| ), we will here 
present the index theorem analysis. Our immediate concern is to count the free parameters 
entering the general soliton solution, with given flux values J2i Qp I d"^ b\r) = ap2nhc{np+ap) 
{p = 1, . . . , M) and the asymptotic behaviors 

const 



oo : l^p(r) 



rpULp 



[ap > 1) . 



(Bl) 



Thanks to Eq. (|3.11| ), the integer Up here coincides with the number of zeros for the field 
\E'p(r). The number of free parameters in question is equal to the number of normalizable 
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zero modes of the small fluctuation equations given in the background of any speciflc soliton 
solution. From Eqs. ( |3.7a| ) and ( p.7b|) , the fluctuation equations are of the following form: 



(B2a) 
(B2b) 



I P' 
[Here, (g^) is being supposed to be a non-singular M x M matrix; but note that our formula 
( [B9| ) is valid for more general (qp)-] Then, to eliminate superflcial zero modes related to the 
freedom of local gauge transformations, we may generalize the real equation ( |B2b| ) to the 
complex equation 



(Vi - 2o-pV2) ?p('^P«i + icTpSai) + 2iap ^ /3pp'^*,5*p/ = . 



(B3) 



Taking the imaginary part of this equation reproduces Eq. ([B2b|) , while the real part can 
be viewed as the gauge condition. 

Equations (|B2a| ) and (p^) are represented by a single matrix equation 



V 



M 



[El qlSa{] + iai [E/ qlSa^2\ 



(B4) 



where 



V 



Di + iaiD2 







-^^1 








Di + i(JMD2 











. 2zfTiAA/^M 


Vi 









• 2iaMpMM'^*M 










(B5) 
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The index of T> is defined as 



Index {T>) = dim(kernel T)) — dim(kernel "D^) 



(B6) 



To calculate this index, it is convenient to consider the quantity 29 



Tr 



VW + M2 



Tr 



M2 



(B7) 



which can be shown to be independent of M^. Naively, one may expect to recover the above 
index in the limit. On the other hand, a straightforward evaluation of I{M'^) in 



the limit cxd gives 



= J2i^p + "p) ■ 



(B8) 



Note that this is not integer-valued in general, while the index defined by Eq. ( |lj6| ) is 
necessarily an integer. This discrepancy is due to the continuum spectrum extending to 
zero, which gives rise to a non-zero contribution to J(M^). Subtracting this contribution 
from the value in Eq. ( p8|) (see Refs. |]20|j30[| ), the correct value for the index reads 



Index {V) = ^(np + dp) 



(B9) 



where dp denotes the largest integer less than a. Further, by the manipulations analogous 
to those described in Refs. pO|,^, it is not difficult to show that dim(kernel "D^) = 0. So the 
kernel of V has the (complex) dimension J2p{i^p + o^p)- Based on this, we now conclude that 
the total number of (real) free parameters in the general solution of the given character is 
equal to 2 J2p{^p + Q^p). Also, in the case of topological soliton solutions which are allowed 
in the presence of non-zero uniform externel magnetic field, it suffices to delete the dp-term 
in this result. 
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FIGURES 

FIG. 1. A closed curve C formed as the positions of two identical particles, r^^-* and r^^-*, are 
exchanged. 

FIG. 2. For parameters lying in the shaded regions, there exists no solution to the self-duality 
equations. 

FIG. 3. Non-Liouville-type solutions in the case of = /?22 = 0. (a)The plot of a solution 
with vorticities n\ = n2 = 0. We have here chosen ^(0) = 0.55 and x(0) = 0.5. (b)The plot of a 
solution with vorticities ni = and n2 = 1. 
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